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A.  EAST  HQ  HAKE  PfiEDICTIOI  THBOfil 

Tbe  desirability  of  being  able  to  predict  when  and  where 
an  earthquake  will  occur  becoa.es  immediately  apparent  in 
light  of  the  recent  devastation  that  occurred  in  Mexico 
City.  The  potential  for  death  if  a  catastrophic  earthquake 
should  cccur  along  the  San  Andreas  fault  is  estimated  to  be 
between  3,000  and  13,000  lost  lives.  The  lower  figure 
applies  for  a  quake  that  hits  at  0230  and  the  higher  figure 
applies  for  a  quake  that  occurs  at  1630  on  a  weekday.  These 
predictions  are  based  on  the  facts  that  wooden,  single 
family  homes  are  much  more  resistant  to  earthquake  struc¬ 
tural  damage  than  are  many  of  the  older,  unreinforced 
masonry  buildings.  In  May,  1983,  30  of  the  40  unreinforced 
masonry  structures  in  Coalinga,  Ca.,  either  partially  or 
totally  collapsed.  In  Los  Angeles,  8,000  similar  buildings 
exist.  The  probability  of  a  magnitude  8.3  earthquake  occur¬ 
ring  along  the  southern  San  Andreas  fault  within  the  next  30 
years  is  estimated  to  be  50  percent.  [ Bef .  1:p.  35] 

The  primary  cause  of  earthquakes  is  the  movement  of  the 
litnosphere,  the  outer  layer  of  the  Earth  consisting  of 
twelve  rigid  plates  riding  on  the  supporting  mantle.  This 
movement  rearranges  continents  and  forms  and  destroys  botn 
mountains  and  oceans.  The  driving  energy  for  this  poten¬ 
tially  large  scale  destructive  mechanism  is  the  convective 
circulation  of  the  mantle.  The  mantle  consists  of  hot, 
solid  rock.  In  time  measured  on  a  geological  scale,  the 
mantle  is  sufficiently  hot  that  it  can  flow  and  be  easily 
deformed.  Hot  magma  from  the  mantle  wells  up  between  sepa¬ 
rating  plates  at  midocean  ridges  forming  new  lithosphere. 
Tne  movement  of  the  cold,  brittle  plates  is  primarily  along 
boundary  faults  and  is  controlled  by  friction.  Earthquakes 


occur  when  these  plates  suddenly  slip  and  elastically 
rebound.  [Ref.  2] 

Figures  1.1  and  1.2  [Ref.  1:p.  38]  show  the  typical 

cycles  of  earthguakes  as  explained  by  the  theory  of  plate 
tectonics.  In  figure  1.1  the  hot  mantle  10-15  km  deep  is 
plastically  deforming  and  stress  is  accumulating  at  asperi¬ 
ties,  or  lock  points,  in  the  cooler  crust,  At  the  moment  of 
rupture  an  asperity  breaks  and  crustal  rupture  spreads 
causing  earthguakes.  The  slippage  is  halted  when  it  encoun¬ 
ters  another  asperity  in  approximately  100  seconds. 
Aftershocks  may  occur  for  the  next  several  months. 
Additional  slippage  at  or  near  the  surface  may  or  may  not 
occur.  In  Figure  1.  2  a  second  major  earthquake  has  occurred 
50  or  100  years  after  the  first.  The  cycle  begins  again 
after  about  200  years  when  the  deep  crust  displacement 
imposes  new  stress  on  the  upper  brittle  crust. 

To  date,  earthquake  prediction  has  been  on  a  long-term 
basis  normally  expressed  as  a  probability  of  quake  occur¬ 
rence  during  the  coming  scores  of  years.  What  has  been 
lacking  is  a  reliable,  short-term  prediction  process. 
Central  to  any  short-term  technique  is  geologic  phenomena 
during  the  period  immediately  preceding  the  main  earthquake. 
Before  the  arrival  of  a  major  shock,  swarms  of  many  small 
earthquakes  are  often  noted  or  sometimes  accelerated  defor¬ 
mation  of  the  crust  is  observed.  However,  neither  one  of 
these  two  predictors  is  reliably  present  in  all  cases. 
[  Ref.  1  :  p.  38  ] 

Systematically  collected  data  on  the  generation  of  elec¬ 
tromagnetic  energy  prior  to,  during,  or  after  an  earthquake 
currently  does  not  exist.  Evidence  of  radio  frequency  (RF) 
emissions  as  a  precursor  to  quake  activity  appears  impul¬ 
sively  suggesting  these  waves  are  generated  over  a  broad 
spectrum  of  frequencies  [Ref.  3],  Citizen  band  and  amateur 
radio  operators  in  the  Hollister,  Ca. ,  area  reported  an 
"increase  in  background  noise  level  preceding  earthquake 
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Figure  1.2  Second  Portion  o£  an  Earthquake  Cycle 


activity"  at  27  MHz  and  at  the  six  and  two  meter  bands. 
This  increase  in  noise  was  on  the  order  of  10  dB  and 
preceded  a  quake  by  12  to  24  hours.  At  one  hour  to  15 
minutes  before  a  quake,  noise  levels  returned  to  normal. 
The  appearance  of  earthquake  lights  [Ref.  4]  and  reception 
of  18  MHz  radio  noise  [Ref.  5]  are  further  examples  of  a 
possible  connection  between  RF  energy  and  earthquakes. 


Various  theories  explaining  the  actual  mechanism  of  BP 
generation  at  the  fault  boundary  include  impact  of 
dissimilar  metals,  shock  induced  polarization  and  fracture 
of  guartz- bearing  rock  [Bef.  6].  Because  of  the  strong 
implication  that  earthquakes  generate  electromagnetic 
energy,  this  thesis  will  present  a  satellite  design  for  a 
space-based  system  as  a  means  of  efficient  quake  detection 
and  rapid  data  dissemination. 

B.  SPACE-BASED  EARTHQUAKE  DETECTION 

The  exploitation  of  outer  space  as  a  tool  in  earthquake 
prediction  based  on  electromagnetic  detection  can  be  divided 
into  two  major  scenarios.  The  first  is  installing  ground- 
based  sensors  that  detect  the  quake  generated  RF  energy  and 
then  relay  data  to  an  orbiting  satellite.  The  satellite 
would  function  as  a  central  collection  point  and  could 
incorporate  features  such  as  data  processing  and  earthquake 
pattern  recognition  with  automatic  alarm  transmission  to  key 
Earth  stations.  The  number  of  ground-based  sensors  to 
provide  adequate  coverage  can  not  be  reliably  estimated 
because  the  strength  of  the  RF  energy  is  not  known  but  is 
being  researched  as  a  portion  of  this  thesis-  As  the 
strength  increases,  the  distance  between  sensors  can 
increase  and  their  number  decrease. 

The  second  scenario  is  direct  detection  of  earthquake 
generated  £F  energy  by  an  orbiting  satellite.  All  the 
previously  mentioned  features  could  be  included.  Direct 
satellite  detection  would  require  a  sufficiently  strong  RF 
signal  to  overcome  the  power  spreading  loss  incurred  from 
the  large  distance  traveled.  The  major  advantages  over  the 
first  scenario  are  worldwide  coverage  with  elimination  of 
the  need  for  ground-based  RF  sensors.  On  an  international 
level,  frequencies  could  be  allocated  to  quake  detection  to 
minimize  possible  man-made  interference  with  the  RF  signal. 
Because  of  the  advantages,  this  second  scenario  will  be  the 
basis  for  the  proposed  satellite  design. 


C.  FOCUS  OF  STUD  I 

A  proposed  design  for  a  space  based  sensor  to  predict 
and  detect  earthquakes  is  presented.  A  free  standing  BF 
noise  measurement  and  recording  system  is  described  to 
research  the  correlation  between  earthquakes  and  increased 
background  electromagnetic  noise  at  30.45  MHz  and  150.75 
MHz.  This  thesis  continues  the  research*  construction  and 
installation  of  equipment  described  in  Reference  6. 
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In  order  to  understand  the  background  information 
that  will  be  used  in  subsequent  portions  of  this  thesis,  it 
is  necessary  to  first  establish  the  coordinate  systea  that 
describes  a  satellite's  position  in  space.  An  epheaeris  is 
a  tabulation  of  the  predicted  or  observed  positions  occupied 
by  the  satellite. 

The  coordinate  systea  to  be  used  is  shown  in  Figure 
2.1  [Bef.  7].  The  large,  flat  rectangular  plane  is  an 
extension  of  the  plane  of  the  equator  and  is  called  the 
equatorial  plane.  Superinposed  on  the  equatorial  plane  is 
the  Cartesian  coordinate  system  indicated  by  T,  T  and  T  with 
the  origin  at  the  Earth's  center.  These  coordinates  do  not 
rotate  but  are  fixed  in  space  with  the  T  axis  pointing  in 
the  direction  of  the  vernal  equinox. 

The  satellite's  motion  about  the  Earth  is  on  the 
circumference  of  the  orbital  plane.  An  elliptical  orbit  is 
shown  but  a  satellite  can  also  have  a  circular  orbit  around 
the  Earth.  The  angular  momentum  or  spin  vector  (h)  for  tne 
satellite  orbit  is  determined  by  the  right-hand  rule  for  a 
south  to  north  equator  crossing  and  is  perpendicular  to  the 
orbital  plane.  Inclination  (i)  is  measured  in  degrees  from 
F  to  the  spin  vector  X.  By  definition,  inclination  can  vary 
from  a  minimum  of  zero  degrees  to  a  maximum  of  180  degrees. 

The  eccentricity  vector  (e)  starts  from  the  origin 
and  points  in  the  direction  of  perigee  or  the  periapsis 
direction.  Perigee  is  the  satellite's  closest  point  of 
approach  to  the  Earth.  The  line  of  nodes  (n)  is  the 
intersection  of  the  equatorial  plane  and  the  satellite's 
orbital  plane.  The  angle  between  the  line  of  nodes  and  the 
eccentricity  vector  is  called  the  arguement  of  perigee  and 
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is  measured  in  the  direction  of  the  satellite  motion. 
Right  ascension  of  the  line  of  nodes  (q)  is  the  angle 
measured  in  the  equatorial  plane  from  the  vernal  equinox 
direction  to  the  line  of  nodes  IT.  The  epoch  angle  (v0  )  is 
the  angle  measured  from  e”  to  the  vector  from  the  center  of 
the  Earth  to  the  satellite  (ra  ).  The  range  of  the  epoch 
angle  is  from  zero  degrees  to  360  degrees. 

2.  Equations  for  aa  Ellipse 

Many  satellites  orbit  the  Earth  in  an  elliptical 
manner  or  use  portions  of  an  elliptical  orbit  during  an 
orbital  transfer  from  one  circular  orbit  to  another  circular 
orbit  of  a  different  radius.  Because  the  equations  of  an 
ellipse  are  an  integral  portion  of  satellite  orbitology,  the 
pertinent  relationships  vill  be  presented.  Figure  2.2  shows 
the  classical  form  of  an  ellipse  with  the  small  circle 
representing  the  Earth  at  one  focus. 

The  semi-major  axis  (a)  is  the  horizontal  length 
measured  from  the  center  of  the  ellipse  to  the  outer  edge. 
The  semi-minor  axis  (b)  is  the  corresponding  vertical 
length.  The  eccentricity  (e)  is  a  measure  of  departure  from 
a  cirlce  and  is  given  by  Equation  2.1 

e  =  SQRT  (1-  (b/a)  2)  (2.1) 

Eccentricity  can  vary  from  a  value  of  zero  indicating  b  =  a 
for  a  circle  to  a  value  of  slightly  less  than  one  if  b 
approaches  zero. 

The  radius  of  apogee  (Ra)  is  the  distance  from  the 
center  of  the  Earth  to  the  point  of  farthest  distance  from 
the  Earth.  The  radius  of  perigee  (Rp)  is  measured  from  the 
center  of  the  Earth  to  the  satellite's  closest  point  of 
approach  to  the  Earth.  The  relationship  between  Ra  and  e 
can  be  expressed  as 

Ra  *  a  (1  *e)  (2.2) 
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tn  to  complete  one  circular  revolution  (Tc)  is  found  by 
dividing  the  distance  traveled  (2rir)  by  the  velocity.  This 
gives 

Te  =  2r*  (2.3) 

v 

=  2tn _  (2.9) 

SQRT (GM/r) 

=  2  t  rV2  (2.10) 

SQRT  (GM) 

Kepler's  third  lav  of  planetary  motion  relates  the 
motion  of  one  planet  to  another.  This  relationship  is  inde¬ 
pendent  of  the  shape  of  their  orbits  and  requires  that  the 
squares  of  the  periods  or  the  time  to  complete  one  revolu¬ 
tion  vary  as  the  cubes  of  the  semi-major  axes.  Then  for  m 
in  an  elliptical  orbit  about  M  with  a  semi-major  axis  of  a, 
the  time  for  one  revolution  (Tc )  is  given  by 

T,  =  2*a3/z  (2.11) 

SQRT  (GM) 

The  elliptical  orbit  velocity  (Ve )  is  given  by 

Ve  =  SQRTfGM  (2/r-1/a)  ]  (2.12) 

2 *  Orbital  Tran  sfers 

It  is  often  desired  to  change  the  orni t  of  a  satel¬ 
lite  to  place  it  either  in  a  more  advantageous  relative 
position  to  other  satellites  or  to  alter  its  relative  posi¬ 
tion  to  the  Eartn  and  allow  it  better  access  to  Earth  gener¬ 
ated  data.  To  perform  an  orbit  modification  the  satellite 
must  expend  energy,  normally  in  the  form  of  burning  fuel 
from  an  on-board  rocket  or  releasing  compressed  gas  tnrough 


Figure  2.3  Gravitational  Force  of  Attraction 


thrusters.  Sockets  are  used  when  a  major  change  is  desired 
because  a  large  force  can  be  guickly  developed  while 
compressed  hydrazine  gas  powered  thrusters  are  normally  used 
for  station  keeping  or  very  minor  orbital  adjustments. 

The  two  major  factors  that  influence  exactly  how  an 
orbit  modification  will  be  accomplished  are  fuel  and  time. 
If  it  is  desired  to  achieve  the  new  orbit  as  rapidly  as 


possible,  then  a  large  amount  of  fuel  will  be  required  to 
provide  the  necessarily  large  force  and  acceleration.  Such 
a  tiae  constraint  nay  exist  in  the  case  of  rescue  aissions 
or  anti-satellite  weapons.  If  time  is  not  the  deciding 
criterion  then  slower,  aore  fuel  efficient  transfers  are 
possible.  This  has  the  aajor  advantage  of  allowing  more  of 
the  satellite's  weight  to  be  in  instr uaentation  since  less 
weight  need  be  allocated  to  fuel.  A  fuel  efficient  transfer 
will  be  used  for  this  thesis  to  naxiaize  operational 
payload. 

Figure  2.4  depicts  a  satellite  at  position  one  in  a 
counter-clockwise  circular  orbit  about  the  Earth.  It  is 
desired  to  discharge  a  payload  from  the  satellite  so  that 
upon  completion  of  the  maneuver,  the  payload  will  be  in 
position  two  and  the  satellite  will  return  to  position  one. 
The  positions  are  separated  by  the  angle  u  measured  in 
radians.  This  same  type  of  maneuver  could  be  used  if  posi¬ 
tion  one  were  to  maneuver  in  crder  to  rendezvous  with  a 
second  satellite  already  at  position  two.  To  accomplish 
either  task,  the  satellite  at  position  one  will  accelerate 
into  the  elliptical  omit  shown.  While  position  two 
continues  on  its  circular  orbit  and  may  pass  the  original 
position  one  several  times,  the  accelerated  satellite  will 
travel  the  elliptical  path  and  rendezvous  will  occur  when 
both  satellite  one  and  two  meet  at  position  one.  Note  that 
the  satellite  at  position  one  can  not  merely  "speed  up"  to 
catch  position  two  because  by  Equation  2.7  they  would  be  at 
different  distances  from  the  Earth.  It  is  necessary  to  know 
how  much  position  one's  velocity  must  be  increased  to  place 
it  into  the  correct  elliptical  orbit  so  rendezvous  will 
occur . 

Since  the  satellites  are  to  meet  at  position  one, 
the  time  for  one  elliptical  orbit  must  be  equal  to  the  tiae 
for  an  integer  number  of  circular  orbits  plus  a  portion  of  a 
circular  orbit  to  account  for  the  initial  angular 


difference.  The  integer  number  of  orbits  (K)  can  be 
selected  to  again  either  minimize  time  or  fuel  expended  with 
more  circular  orbits  indicating  less  fuel  required.  That  is 

KTC  +  Tc(1-u/2ir)  =  Te  (2.13) 

Substituting  Equations  2.10  and  2.11  into  Equation  2.13 
gives 

K21rRe3/2  *  2irEc  3/2  (1-U/2^)  =  2  Td  3/2  (2.14) 

SQfiT(GM)  SQRT(Gd)  SQRT  (GM) 

The  circular  orbit  radius  is  now  being  represented  by  Re 
rather  than  just  r  to  assist  in  equation  development.  The 
common  term  of  2ir/S2RT(GM)  can  be  cancelled  which  leaves 

KRC  3/2  +  fie  3/2  (  |-u/2^>  =  a3/2  (2.15) 

Rearranging  and  solving  for  a  gives 

a  =  Rc  [K  +  l-u/2  ir]2/3  (2.1b) 

Eguation  2.16  gives  an  expression  for  a  which  is  tne 
semi-major  axis  of  the  required  elliptical  orbit.  An 
expression  is  needed  that  relates  the  eccentricity  of  the 
transfer  ellipse  to  K.  This  can  be  done  by  using  Equation 
2.3  and  realizing  that  at  perigee  Rc  =  Sp.  taxing  the 
substitution  yields 

a  =  a  (1-e)[  K+l-u/2  tt]2/3  (2.17) 

Cancelling  the  a  and  rearranging  gives 

1  =  [  K  +  1-u/2ir]2/3  -  e[  KH-u/2  tt]2/3  (2.18) 


Solving  for  e  gives  Equation  2.19 


e  *  [  K+1-u/2ir]2/3  -  1  (2.19) 

[  K  +  l-u/2  tt]2 /3 

Equation  2.19  now  relates  the  necessary  eccentricity 
of  tbe  transfer  ellipse  to  both  the  number  of  circular 
orbits  of  position  two  and  to  their  initial  angular  separa¬ 
tion.  Equation  2.7  gives  the  circular  orbit  velocity  and 
Equation  2.12  describes  the  elliptical  orbit  velocity.  The 
velocity  of  the  elliptical  orbit  at  Rp  must  be  known  to  find 
the  difference  between  the  two  velocities.  Evaluating 
Equation  2.12  at  r  =  Rp  gives 

Vp  =  SQET[GM(2/Rp- 1/a)  ]  (2.20) 

Vp  is  the  velocity  at  perigee  for  satellite  one.  Solving 
Equation  2.3  for  a  gives 

(2.21) 

Substituting  Equation  2.21  into  Equation  2.20  and  simpli¬ 
fying  yields 

Vp  =  SQRT[  (GH/Rp)  (He)  ]  (2.22) 

The  required  change  in  velocity  can  now  be  computed 
by  subtracting  Ve  from  Vp. 

Vp  -  vc  =  SQRI[  (GH/Sp)  (He)]  -  SQRT[  GM/Re  ]  (2.23) 

At  perigee  fle  =  Rp;  Vpe  will  represent  the  difference  in 
velocity. 

Vpc  =  SQfeT  (GM/R<)  [  SQRT  (He)  -  1]  (2.24) 


The  square  root 


of  (1  ♦  e)  can  be  approximated  by 


1  *'.5e  ♦  .  .  .  .  Also  note  that  SQBT(GH/He) 
Equation  2.24  can  be  written  as 


Thus 


-  ve  - 


Vpc  ~  ve  [  1  ♦  .  5e  -  1] 


(2.25) 


■  Yce  (2.26) 

2 

=  Vc[  (K+1-U/2*-)  2/3  -  1  ]  (2.  27) 

2[K*1-U/2ir]2/3 

Equation  2.27  is  a  key  formula  and  can  be  used  to 
determine  how  much  additional  velocity  is  required  to  ensure 
a  proper  rendezvous  transfer  orbit.  It  has  the  flexibility 
to  vary  the  number  of  circular  orbits  which  is  equivalent  to 
being  able  to  predetermine  the  rendezvous  time. 

The  use  of  this  equation  can  be  extended  to 

deploying  numerous,  equally  spaced  satellites  from  a  single 
large  parent  vehicle.  The  parent  vehicle  would  be  in  the 
elliptical  orbit  carrying  the  satellites  to  be  deployed. 

Each  time  the  "rendezvous  point"  is  achieved  a  satellite 
would  be  ejected.  Two  satellites  would  not  be  placed  in  the 

same  position  as  the  new  rendezvous  point  is  always  u 

radians  ahead  of  the  current  one.  Because  Equation  2.27  has 
u  as  a  variable  parameter,  any  angular  positioning  is 
possible  and  numerous  satellites  may  be  deployed  limited 
only  by  the  total  weight  that  can  be  launched  from  Earth. 

3.  Orbital  Transfer  Fuel  Requirements 

The  next  parameter  that  must  be  calculated  is  how 
much  fuel  is  required  to  inject  the  parent  satellite  from 
the  circular  orbit  into  the  elliptical  orbit.  Comparison  of 
Equations  2.7  and  2.22  show  that  Vp  has  a  larger  magnitude 
than  Vc  even  though  both  are  at  a  distance  Sp  from  the 
center  of  the  Earth.  Thus,  a  transfer  from  a  circular  to 


elliptical  orbit  requires  the  satellite  to  accelerate  in  the 
direction  of  travel.  Similarly,  when  the  parent  satellite 
ejects  a  payload,  that  payload  must  burn  fuel  and  slow  down 
to  Vc  or  it  would  continue  on  the  parent's  elliptical  path. 
The  amount  of  fuel  required  for  either  maneuver  is  critical 
since  not  only  the  fuel  but  also  its  container  subtract  from 
the  weight  that  could  be  used  for  instrumentation. 

Rochet  fuel  energy  is  classified  by  its  specific 
impulse  <Isp)  and  is  a  measure  of  the  thrust  developed  as  a 
function  of  fuel  mass  flow  rate.  The  mass  flow  rate  can  be 
in  terms  of  gallons  per  minute  for  a  liquid  fuel  or  pounds 
per  second  for  a  burning  solid  fuel.  The  accepted  English 
units  for  I$p  are  expressed  as  seconds.  This  comes  from 


=  pounds  of  thrust 

pounds  of  fuel/second 


=  seconds 


(2.  28) 


Typical  values  of  Isp  range  from  275-325  seconds.  Knowing 
the  fuel  Isp,  the  force  of  acceleration  or  thrust  (Ff)  can 
be  found  by 


F,  =  Isp  (dfci/dt) 


(2.29) 


H  is  the  propellent  weight  flow  rate. 

The  rocket  thrust  is  generated  from  burning  fuel 
that  exits  the  rocket  at  a  velocity  Vx  with  units  of  feet 
pec  second.  T ne  fuel  Isp  and  exhaust  velocity  can  oe 
related  by 


V*  =  gi*. 


(2.30) 


where  g  is  the  acceleration  constant  equal  to  32.2  ft/sec2. 
Tne  thrust  can  also  te  related  to  the  exhaust  velocity  and 
fuel  flow  rate  by 


F,  =  -  (dm/dt)  V. 


(2.31) 


The  minus  sign  is  necessary  because  the  exhaust  is  traveling 
in  one  direction  and  the  resultant  force  is  in  the  opposite 
direction.  By  Newton's  laws,  the  rocket  must  react  to  the 
fuel  thrust  by 


F  =  M$  a  (2.32) 

hs  is  the  mass  of  the  rocket  being  accelerated  and  a  is  the 

resultant  acceleration.  The  acceleration  can  be  written  in 
terms  of  the  rocket  velocity  (Vr)  as 

a  *  dVr/dt  (2.33) 

Substituting  Eguation  2.33  into  Equation  2.32  gives 

F  =  Ms  (dvr/dt)  (2.34) 

Ft  must  equal  F.  Setting  Equation  2.34  equal  to  Equation 
2.31  gives 


M,  (dvr  /dt)  =  -  (dm/dt)  V, 


(2.35) 


Simplifying  and  rearranging  terms  gives 


dVr  =  -Vx  (da/Ms ) 


(2.36) 


The  left  side  of  this  equation  can  now  be  integrated  from 
the  initial  rocket  velocity  ( V;  )  to  the  final  rocket 
velocity  (Vj  )  .  The  right  side  can  oe  integrated  from  the 
initial  rocket  weight  (flj )  to  the  final  rocket  weight  ( ) . 


W, 


fr  =  -Vx 

i  W; 


(2.37) 


Equation  2.38  shows  the  evaluation  result 


vf  -Vj  =  -Vxln(Hf /W;)  <2.38) 

Dividing  by  -Vx  and  using  both  sides  as  a  power  of  fche  expo¬ 
nential  e  gives 

exp£-(Vf-V;  )/Vx  ]  =  tif/W;  (2.39) 


V,  =  *i  exp[-(V, -V;  )/Vx  ]  (2.40) 

The  initial  weight  oust  be  the  sun  of  the  final  weight  plus 
the  weight  of  any  fuel  or  propellent  that  was  burned  (Wp)  to 
produce  the  change  in  velocity,  or 

'4,  =  Wj  -lip  (2.  41) 

Substituting  Equation  2.41  into  Equation  2.40  gives 

H;  -Wp  *  Hj  exp(  -  (Vf-Vj  )  /Vx  ]  (2.42) 

Solving  for  Hp  gives 

Wp  =  Hj  [  1  -  exp-(Vf-Vj  )/Vx  ]  (2.43) 

Replacing  Vx  with  glsp  from  Equation  2.30  gives 

Wp  =  Wj  [  1  -  exp-(Vf  -Vi  ) /(glsp)  ]  (2.44) 

The  quantity  Vj-Vj  is  the  same  change  in  velocity  required 
by  Equation  2.27  for  orbital  transfer. 


VpC  =  Vf-Vi 


(2.45) 


Equation  2.44  then  becomes 


»P  -  Wit  1  -  exp-<ypc/gl#p)  ] 


(2.46) 


Equation  2.27  can  be  used  to  deteraine  tbe  amount  of 
velocity  change  that  is  needed  for  an  orbital  transfer  and 
Equation  2.46  is  used  to  deteraine  the  percentage  of  initial 
satellite  weight  that  must  be  reserved  for  the  fuel  to 
accomplish  the  same  transfer. 
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A.  GEHEEAL  DESIGM 

The  satellite  will  be  placed  in  a  circular  low  Earth 
orbit  (LEO)  of  322  lea.  An  orbit  close  to  the  Earth  will  be 
used  because  of  the  anticipated  low  power  level  of  the 
earthquake  electromagnetic  signal.  The  design  of  the  satel¬ 
lite  will  be  based  on  using  the  space  shuttle  for  achieving 
the  initial  circular  orbit.  Because  of  the  large  cargo  bay 
on  the  shuttle  (60  ft  X  15  ft)  ,  a  group  of  four  satellites 
arranged  in  a  circle  and  mounted  on  a  central  support  bus 
will  comprise  the  launch  package.  When  deployed,  the 
satellites  will  be  egually  spaced  by  90  degrees. 

Primary  electrical  power  will  be  provided  by  solar  cells 
with  batteries  supplying  power  during  the  time  the  satellite 
is  in  the  Earth’s  shadow.  Spin  stabilization  of  the  satel¬ 
lite  body  with  hydrazine  thrusters  for  station  keeping  will 
be  used  to  maintain  a  de-spun  antenna  pointing  at  the  Earth 
with  attitude  sensing  to  be  done  by  Earth  limb  sensors. 
Each  satellite  will  weigh  approximately  500  pounds  and  have 
a  cylindrically  shaped  main  body- 

B.  A  ST  EMI  A  QESIGM 

The  satellite  will  be  designed  to  detect  electromagnetic 
energy  in  the  30  MHz  region.  This  frequency  is  used  because 
it  is  designated  for  astronomical  observations  and  manmade 
interference  should  be  at  a  minimum.  Additionally  30  MHz  is 
well  above  the  cutoff  frequency  (Fp)  for  propagation  through 
the  ionosphere  given  by 

Fp  =  9SQfiI(M.)  (3.1) 

where  N*  is  the  ionospheric  plasma  electron  density  with 
units  of  electrons/meter3.  A  nominal  value  for  Fp  is 


10  HHz.  Any  frequency  higher  than  Fp  will  successfully 
penetrate  into  outer  space  with  lower  frequencies  being 
reflected  bach  to  Earth. 

A  Yagi  antenna  will  be  used  as  the  main  sensor  because 
of  the  size  and  space  constraints  imposed  by  a  satellite 
design.  Refer ence  8  will  be  used  as  a  basis  for  the  design 
of  the  antenna.  Table  1  [Ref.  8s  p.  7]  is  used  once  the 
initial  design  criteria  is  established.  This  initial  data 
consists  of  the  frequency  or  operation  (30  MHz) ,  antenna 
gain  desired  (9.2  dB )  ,  d/A  ratio  where  d  is  the  diameter  of 
the  parasitic  elements  (2  cm),  and  the  D/a  ratio  where  D  is 
the  diameter  of  the  support  boom  (4  cm)  .  The  wavelength  A 
is  given  by 

A  =  c  (3.2) 

f 

Tne  speed  of  light  (c)  is  a  constant  equal  to 
3  X  10®  meters/second  and  f  is  the  operating  frequency. 
Thus  a  frequency  of  30  MHz  gives  a  wavelength  of  10  meters. 
Following  the  steps  in  Reference  8,  Figure  3.1  [Ref.  8:p. 
20]  gives  the  length  of  the  first  and  third  director  as 
4.44  meters,  4.5  meters  for  the  length  of  the  second 
director  and  4.9  meters  for  the  length  of  the  reflector. 
The  element  spacing  is  2  meters.  Deployment  of  the  antenna 
in  outer  space  can  be  achieved  by  using  a  bistem  which  is  a 
deployable  boom  rolled  up  on  a  drum  in  its  flattened  config¬ 
uration.  Because  of  its  small  size  when  stowed,  total 
storage  space  required  will  be  approximately  one-half  a 
cubic  foot.  Since  the  antenna  will  be  operating  in  a 
weightless  environment,  strong  structural  members  are  not 
required.  The  total  weight  of  the  antenna  would  not  exceed 
30  pounds  since  the  material  of  construction  could  be  a 
light-weight  aluminum  based  alloy.  The  gain  of  9.2  dB  is 
referenced  to  a  dipole  antenna.  For  reference  to  an 
isotropic  radiator,  2.16  dB  of  gain  must  be  added. 
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C.  FUEL  REQUIRES  BITS 
1 .  Satellite  Fu el 

The  on-station  weight  of  the  satellite  is  590  pounds 
and  this  weight  must  be  decelerated  from  the  parent1 s  ellip¬ 
tical  velocity.  Equation  2.27  will  be  used  to  determine  the 
necessary  velocity  change  and  Eguation  2.7  will  be  used  to 
find  Vc .  GM  evaluates  as  a  constant  equal  to 
3.99  X  10l*  m3/sec2  and  r  is  the  sum  of  the  radius  of  the 
Earth  (6370  Km)  and  the  satellite  altitude  above  the 
Earth  (32 2  Km). 


=  SQRT  (3. 99  X  10^/6692  X  10^) 


(3.3) 


=  7,722  m/sec 


=  25,335  ft/sec 


(3.4) 


(3.5) 


Equation  2.27  can  now  be  evaluated  if  K  is  set  equal  to  one 
and  u  is  set  equal  to  ir/2  radians. 


Vpc  =  25, 335[  (2-.2S)2/3  -  1] 

2[  (2-.  25)  a/3  -  i  ] 


=3,945  ft/sec 


(3.6) 


(3.7) 


Equation  2.46  can  be  used  to  find  the  fuel  fraction  given  oy 


WD/Wi 


1  -  exp-r  3945/(32)  (300)  ] 


(3.8) 


CllniNT  DIAMETER 


The  initial  weight  is  composed  of  the  cn- orbit 
satellite  weight  of  500  pounds  (H, )  ,  and  the  weight  of 
the  case  (We)  that  contains  the  fuel.  That  is 


Hj  = 


(3.10) 


For  an  efficiently  designed  rocket  motor,  the  majority  of 
the  weight  is  due  to  fuel  with  the  remaining  weight  due  to 
the  case.  This  ratio  (f)  is  given  by 


£  = 


(3.11) 


Wp  ♦Wc 


A  typical  value  of  f  is  0.93.  Using  this  definition  of  f. 
Equation  3.10  can  be  rewritten  as 


W i  =  »,  ♦Hp/f 


(3.  12) 


Substituting  Wp  =  ,337Wj  from  Equation  3.9  and  0.93  for  f 
into  Equation  3.  12  gives 


Hi  =  W,  ♦  .  337 Kj  /.93 


(3.  13) 


Rearranging  terms  to  solve  for  Wj  and  substituting 
500  pounds  for  Ws  gives 


»i  =  500 


1  -  .337/. 93 


(3.  14) 


=  784  lb 


(3.  15) 


This  can  be  used  to  find  Hp  from  Equation  3.9 


WB  =  (784)  (.337) 


(3.  16) 


Wp  =  264  lb 


(3.17) 


The  Height  of  the  case  can  be  found  by  rearranging 
Eguat ion  3.11 


Bc  =  Hp(1-f) 
f 


(3.  18) 


=  264  (1-.93) 


(3.19) 


=  20  lb 


(3.20) 


2  .  Parent  Satellite  Fuel 

The  amount  of  fuel  required  by  the  parent  to  inject 
itself  into  the  elliptical  orbit  can  now  be  calculated  in  a 
similar  manner.  The  central  support  bus  weight  will  be  set 
at  500  pounds.  This  relatively  large  weight  is  required 
because  the  bus  must  be  strong  enough  to  rigidly  hold  the 
four  satellites  during  the  shuttle  launch.  Thus,  enough 
fuel  must  be  allowed  to  accelerate  the  bus  and  three  satel¬ 
lites  into  an  elliptical  orbit  since  one  satellite  is 
deployed  prior  to  the  parent's  acceleration.  The  net  weight 
to  be  boosted  is  then  2,852  pounds.  Similar  to  Equation 
3.14,  the  initial  total  weight  is 


!*i  *  2852 _ 

1  -  .337/. 93 


(3.21) 


=  4,473  lb 


(3.  22) 


The  fuel  weight  and  case  weight  for  the  parent  can  also  be 
found  in  the  same  manner  as  for  the  individual  satellite. 


w  =  (4473)  (.33  7) 


(3.  23) 
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Hp  =  1,507  lb  (3. 24) 


Hc  =  1507(1-. 93)  (3.25) 


=  114  lb  (3.26) 


D.  ELECTRICAL  POSER 

The  primary  source  of  electrical  power  will  be  provided 
by  converting  light  energy  into  electrical  energy  by  use  of 
solar  cells.  Because  the  available  power  from  an  individual 
cell  is  small  (on  the  order  of  0.06  watts)  ,  a  large  number 
of  cells  must  be  used  that  collects  solar  energy  over  a 
large  incident  area.  In  order  to  provide  the  desired  total 
power  at  useful  voltages,  solar  cell  configurations  are 
normally  used  in  a  series-parallel  connection. 

The  satellite  will  not  be  in  sunlight  during  the 
Earth-Sun  eclipse  periods  of  the  orbit.  During  these  times, 
rechargeable  batteries  will  provide  the  electrical  power. 
Additionaly,  the  batteries  can  provide  bach-up  power  during 
transients  that  exceed  the  solar  panel  capacity. 

In  order  to  provide  a  useful  power  system,  these  two 
power  sources  must  be  integrated  with  a  power  control  and 
power  conversion  system.  Functions  that  must  be  accom¬ 
plished  include  battery  charging  during  sunlight  periods, 
bus  voltage  monitoring  and  control  and  converting  battery 
power  to  electrical  power  when  required. 

Figure  3.2  [Ref.  9:p.  1-7.5]  shows  a  block  diagram  of  a 
satellite  power  system.  While  only  one  of  each  component  is 
shown,  normal  installations  consist  of  multiple  components 
for  reliability  reasons.  The  power  from  the  main  solar 
panel  supplies  the  spacecraft  systems  through  the  primary 
power  distribution  bus.  The  output  voltage  of  the  panel  can 


▼ary  as  a  function  of  load,  temperature  and  other  factors. 
Bus  voltage  limiters  prevent  the  voltage  swings  from 
exceeding  spacecraft  design  ranges  by  acting  as  variable 
shunt  loads  during  a  rising  voltage  condition.  The  battery 
discharge  regulator  turns  on  when  a  pre-determined  bus  low 
voltage  is  sensed  and  provides  a  regulated  output.  The 
battery  terminal  voltage  is  higher  than  the  required  bus 
voltage  and  in  order  to  charge  the  battery  a  still  higher 
level  of  voltage  is  needed.  To  provide  this  higher  voltage 
level,  the  additional  solar  panel  is  series  connected  to  tne 
main  panel.  The  size  of  the  battery  charging  panel  is 
designed  so  that  the  charge  control  functions  provided  in 
the  controller  need  only  be  an  on-off  and  rate  change 
switch.  The  voltage  and  current  sensors  provide  means  to 
monitor  the  operational  state  of  the  power  system  and  for 
load  managing  during  the  orbit. 

1.  Solar  Panel  Design 

As  shown  in  Figure  3.3  £Bef.  9:p.  1-7.11]  individual 
solar  cells  are  arranged  in  a  series-parallel  array  to 
provide  the  required  power.  The  three  strings  are  in 
parallel  to  minimize  power  loss  should  a  single  cell  fail. 
These  arrays  are  mounted  on  the  satellite  with  the  series 
string  length  parallel  to  the  satellite  spin  axis.  All  the 
arrays  are  then  connected  in  parallel  to  produce  the  desired 
power.  The  surface  of  the  solar  cell  is  considered  to  be  a 
Lambertian  surface  and  the  power  produced  by  the  cell 
decreases  by  a  cosine  function  as  the  sun's  incident  angle 
changes  from  normal  as  shown  in  Figure  3.4  [Bef.  9:p. 
1-7.  11], 

Figure  3.5  £Bef.  9:p.  2-4.17]  shows  how  the  sun's 
incident  angle  varies  during  one  year.  When  the  solar  cells 
are  in  the  sunlight,  only  those  positioned  in  the  center  are 
illuminated  by  normal  rays.  Going  left  or  right  from  center 
causes  the  cell  output  to  drop  by  the  cosine  function.  The 
net  result  is  an  output  from  the  half-cylinder  that  is  equal 
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Figure  3.3  Solar  Cell  Array 

to  the  output  from  a  flat  plate  whose  dimensions  are  the 
diameter  and  height  of  the  cylinder.  Additionally,  the  sun 
line  varying  23.5  degrees  as  shown  in  Figure  3.5  causes  a 
further  decrease  in  solar  panel  output  power. 

Solar  cell  performance  also  degrades  with  time  from 
exposure  to  electron  and  proton  fields  trapped  in  the 
Earth’s  magnetic  field  and  solar  flares.  The  extent  of 
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Figure  3.4  Cylindrical  Panel  Installation 


degradation  is  shown  is  Figure  3.6  [Ref.  9:  p.  1-7.13]. 
Danage  is  a  result  of  charged  particles  creating  defects  in 
the  n  and  p  layers  of  the  cell.  The  net  result  is  shown  in 
Figure  3.7.  For  any  set  of  four  sequential  peaks,  the  maxi- 
nuns  correspond  to  the  autumnal  and  vernal  equinox  and  the 
minimums  correspond  to  suamer  and  winter  solstice  when  sun 
incidence  is  offset  by  23.5  degrees. 


Figure  3.5  Solar  Incident  ingle 


The  satellite  that  is  being  proposed  will  have  a 
diameter  of  5.7  feet  and  a  height  of  10  feet.  This  gives  a 
projected  area  of  57  square  feet.  At  beginning-of-lif e 
(BOL)  solar  cells  produce  approximately  11  watts  of  power 
per  sguare  foot  of  area.  Thus  the  57  square  feet  will  yield 
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electrical  power  budget  Bust  be  based  on  this  EOL  value. 
The  power-to-mass  ratio  at  EOL  for  solar  cells  is  approxi¬ 
mately  20  vatts/kg  for  normal  incidence  [Bef.  10:p.  73]. 
This  ratio  must  be  decreased  because  the  satellite  does  not 
have  flat  plate  collectors.  The  amount  of  decrease  will  be 
the  ratio  of  the  area  of  a  flat  plate  collector  to  the  area 
of  the  illuminated  half  cylinder.  The  area  of  an  eguivalent 
flat  collector  is  57  sguare  feet.  The  surface  area  of  the 
half  cylinder  (Ac )  is 


Ac  =  (-5)  (2n-rh) 


-  (it)  (5.7/2)  (10) 


=  89.5  ft* 


(3.28) 


(3.29) 


(3.30) 


The  power- to- mass  ratio  of  the  cylinder  (Pe )  is  then 


Pc  =  (20)  (57/89.5) 


=  12.  73  watts/kg 


(3.31) 


(3.  32) 


The  mass  of  the  solar  cells  (Mc)  is  then 


Mc  =  200/12.73 


=  15.7  kg 


(3.33) 


(3.  34) 


This  is  the  mass  of  only  those  cells  that  are  receiving 


sunlight.  The  total  mass  is  twice  as  ouch.  Converting  kg 
to  pounds  will  give  a  total  weight  of  the  solar  cells  as 
69  pounds. 

2 .  Batter v  Reguireo ents 

Batteries  that  are  used  in  a  weightless  environment 
can  not  be  of  the  "flooded"  electrolyte  type  such  as  lead- 
acid  cells;  Rather,  the  electrolyte  oust  be  absorbed  in  the 
separator  material  to  assure  even  distribution.  The  two 
primary  types  of  batteries  that  meet  this  condition  are  the 
sealed  silver-zinc  and  the  nickel-cadmium. 

Table  II  £Ref.  9:p.  1-7.17]  shows  the  major  charac¬ 
teristics  of  concern  for  both  battery  types.  Nickel-cadmium 
batteries  will  ne  used  because  of  their  longer  life  expec¬ 
tancy  and  tolerance  to  overcharge  allowing  a  simpler  charge 
control  circuit.  The  major  disadvantage  will  be  an  increase 
in  weight  due  to  the  lower  energy  density. 

The  satellite  will  be  spending  50  percent  of  its 
time  in  the  Earth's  shadow  during  which  tne  battery  must 
supply  the  electrical  power.  During  the  time  sunlight  is 
available,  the  battery  must  he  recharged.  Since  equal 
amounts  of  time  are  spent  in  light  and  shade,  half  of  the 
solar  panels  must  be  dedicated  to  battery  charging  and  the 
remaining  half  to  satellite  bus  power.  Ihe  time  the  satel¬ 
lite  is  in  darkness  (Td)  can  be  found  from  Equation  2.10  to 
be 


Td  =  .5TC 


(3.35) 


=  (7r)  (b692  X  10  3)  3,/2/sqrt  (3.99  X  10**)  (3.36) 


T 


45.4  min 


<3.  371 


ENERGY 


The  satellite  system  vill  be  drawing  200  watts  of  power 
during  this  time  and  the  watt-hours  required  from  the 
battery  (?b)  will  be 


Pb  =  (200)  (4  5.4/60) 


(3.38) 


»  151.3  watt-hours  (3.39) 

Figure  3.8  fBef.  9:p.  1-7.21]  shows  taat  as  the 
depth  of  discharge  (DOD)  decreases  the  lifetime  of  the 
battery  increases.  In  order  to  have  a  long  lifetime,  a  high 
capacity  battery  system  is  desirable  but  the  penalty  is 
additional  weignt.  From  Reference  10,  the  battery  system  on 
the  Intelsat  V  can  supply  34  amp-hours  at  27  volts  and  has  a 
mass  cf  32  kg  [Ref.  10:p.  73].  This  gives  a  battery  energy 
density  of  13  watts/pound  which  is  slighty  better  than  shown 
in  Table  II.  To  provide  a  low  DOD  but  not  require  excessive 
weight,  the  battery  system  for  the  proposed  satellite  will 
be  capable  of  supplying  1200  watt-hours.  This  gives  a  DOD 
of 


DOD  =  15  1.3/1200 


(3.40) 


=  12.o  percent  (3.41) 

From  Figure  3.9  this  gives  approximately  23,000  cycles  for 
the  battery.  The  time  for  one  cycle  is  90.8  minutes.  There 
are  5,789  cycles  completed  in  one  year.  Thus  the  battery 
has  a  lifetime  of  approximately  four  years  and  will  weigh 
about  92  pounds. 


E.  SIS1EM  HEIGHT  AH D  POSES  BUDGET 


The  satellite  has  an  on-orbit  weight  of  500  pounds. 
Subtracting  30  pounds  for  the  antenna,  69  pounds  for  the 
solar  cells  and  the  battery  weight  of  92  pounds  leaves 
309  pounds  for  the  remaining  electronics,  structure  and 
control  systems.  Using  the  weight  and  power  ratios  in 
Reference  10  gives  the  results  shown  in  Table  III. 


TABLE  III 

HEIGHT  AHD  POHER  BUDGET 


System 

Height  (lb) 

Power (wat 

Antenna 

30 

— 

Solar 

Cells 

69 

— —  — — * 

Battery 

92 

Telemetry  , 
Command 

16 

7.9 

Attitude 

43 

9.6 

Control 

Propulsion 

24 

1.0 

Structure 

86 

— 

Electrical 

Integration 

25 

20.3 

Signal 
Processin  g 

115 

161.2 

Of  the  161.2  watts  available  for  signal  processing 
approximately  30  percent  would  be  available  for  use  t 
transmit  data  back  to  central  Earth  collection  stations 
The  signal  strength  at  the  surface  of  the  Earth  can  b 
calculated  by 

Pr  =  ?t  (GtGr) 

(4*)  2C2 


(3.  42) 


Pt  is  the  transmitted  power,  Gt  and  Gr  are  the  gains  of  the 
transmitting  and  receiving  antennas  and  r  is  the  distance 
between  the  satellite  and  Earth.  Using  the  designed  antenna 
as  the  transmitting  antenna,  the  gain  is  11.36  dB  with 
reference  to  an  isotropic  radiator.  An  Earth  based  antenna 
does  not  have  the  same  weight  constraints  and  could  be 
designed  for  a  gain  of  150.  Osing  these  values  ,  tne 
received  signal  power  would  be  0.63  microwatts. 

F.  STABILITY 

The  satellite  is  designed  to  spin  about  its  longitudinal 
axis.  If  it  were  to  tumble  about  its  transverse  diameter, 
it  would  be  totally  inoperative.  An  analysis  of  the  moment 
of  inertia  about  both  axes  is  required  to  ensure  proper 
placement  of  the  internal  components.  From  Reference  11  the 
moment  of  inertia  for  a  hollow  circular  cylinder  of  length  h 
about  the  longitudinal  axis  is 


Ig  =  (m/2)  (R  j  2*R02) 


(3.43) 


Rj  and  R0  are  the  inner  and  outer  radii  and  m  is  the  mass  or 
the  cylinder.  The  moment  of  inertia  about  the  transverse 
diameter  is 


It  =  m[R;2+R02  f  h2J 


(3-  44) 


For  stability,  Ig  must  be  greater  than  lt.  The  longitudinal 
moment  of  inertia  due  to  the  solar  cells  with  Rj=R0  yields 
560.4  ft2- lb.  The  corresponding  transverse  moment  or 
inertia  is  655.2  ft2-lb.  In  order  to  increase  Ig  more  than 
It,  weight  must  be  distributed  in  a  disx-like  manner  at  the 
center  of  the  satellite.  The  weight  of  the  antenna,  solar 
cells  and  structure  can  not  be  centralized.  The  remaining 
weight  is  315  pounds.  This  weight  with  Rj  equal  to  zero 
adds  1279.3  ft2-lb  in  the  longitudinal  direction  for  a  total 


of  t 839.7  ft2-lb.  It  is  now  necessary  to  determine  the 
maximum  height  (d)  of  the  disk.  The  total  transverse  moment 
of  inertia  must  be  less  than  1839.7  ft*-lb.  That  is 


IV.  DESCRIP TIP B  g£  BESE1BCH 


A.  DESIGH  OVERVIEW 

The  earthquake  monitoring  system  is  installed  three 
miles  east  of  San  Juan  Batista,  ca.,  and  seven  miles  south 
of  Hollister.  It  is  located  in  the  Hollister  Hills  State 
Vehicular  Recreation  Area  (Upper  Ranch)  south  of  the  MX 
track  (see  Figures  4.1  and  4.2).  The  site  was  chosen 
because  it  is  both  historically  and  recently  an  active 
earthguake  area  and  is  far  enough  from  industrial  areas  to 
minimize  manmade  noise  signals  yet  near  possible  earthquake 
sources  to  generate  electromagnetic  energy.  Exact  location 
is  1000  feet  from  the  San  Andreas  fault  at  121°  23.5’  w  and 
360  45.5*  N. 

Two  frequencies  (38.45  MHz  and  1  50.75  MHz)  are  being 
monitored.  These  frequencies,  each  with  a  50  KHz  bandwidth, 
have  been  assigned  to  this  project  by  the  Army  Frequency 
Coordinator  for  the  Western  United  States  for  a  renewable 
one  year  period.  The  antenna  complex  consists  of  two 
20  foot  towers  and  mounted  antennas  as  shown  in  Appendix  D. 
Each  tower  has  five  mounted  antennas.  Eight  Yagi  antennas 
sense  horizontal  and  vertical  polarization  at  each  frequency 
both  up  and  down  the  fault.  Additionally,  an  omni  antenna 
at  each  frequency  is  vertically  mounted. 

Each  antenna  signal  is  individually  cabled  to  a  modular 
AM  receiver.  Output  from  the  receiver  is  sent  to  an  inter¬ 
face  amplifier  where  it  is  joined  with  a  system  timing 
signal.  The  composite  signal  is  then  sent  to  a  chart 
recorder  and  permanently  traced  onto  pressure  sensitive 
paper.  Figure  4.3  shows  a  block  diagram  for  each  channel 
with  a  single  timer  signal  being  common  to  all  ten  channels. 
All  components  are  powered  by  12  V  DC.  Collected  data  will 
be  correlated  to  earthquake  data  provided  by  the  Geological 
Survey  Office,  Menlo  Park,  Ca. 
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Fig  are  4.1  Site  Location 


B.  SYSTEM  COBPOB ENTS 

1 .  fiecei vers 

The  antenna  signal  is  coupled  into  the  RF  amplifier 
by  a  50  ohm  RCA  jack.  The  oscillator  frequency  is  deter¬ 
mined  by  the  crystal  Y1.  For  38.45  MHz ,  the  crystal 
frequency  is  49.15  MHz  and  for  150.75  MHz  the  crystal 
frequency  is  46.68333  MHz.  Figure  4.4  shows  the  schematic 
diagram  for  the  150.75  MHz  receiver.  The  38.75  MHz  receiver 
does  not  have  the  Q4  multiplier  circuit. 


Figure  4.2  Site  Location 


The  10.7  MHz  IF  signal  is  sent  to  the  IF-20  IF/Audio 
assembly.  The  signal  is  mixed  with  the  10.245  MHz  oscil¬ 
lator  frequency  and  down-converted  to  455  KHz.  This  signal 
is  then  amplified  and  detected.  The  MS"  meter  output  from 
an  integral  log  detector  is  used  as  the  input  to  the  inter¬ 
face  amplifier.  Figure  4. 5  is  a  schematic  diagram  for  the 
IF/detector  circuit.  Appendix  B  contains  measured  receiver 
sensitivity  curves. 

2.  Timer 

In  order  to  determine  the  time  of  data  reception, 
the  timer  generates  a  voltage  signal  every  15  minutes, 
hourly  and  every  24  hours.  The  duration  of  this  signal  was 
modified  from  the  original  design  to  new  values  of 
14  seconds  for  the  15  minute  pulse,  42  seconds  for  the 
hourly  pulse  and  98  seconds  for  the  24  hour  pulse.  Accurate 
time  is  reguired  to  allow  comparison  with  earthquake  data. 
A  21  stage  counter  timer  (RCA  CD4045A)  is  crystal  controlled 
and  establishes  a  one  second  time  reference  for  the  entire 
system.  Synchronous  down  counters  (RCA  CD4010  3B)  then  use 
this  reference  to  generate  signals  at  the  appropriate  times. 
The  reset  button  is  used  to  initialize  or  reset  the  zero 
time  mark.  A  clock-runniag  LED  is  used  as  a  visual  indica¬ 
tion  of  proper  operation  and  a  pen-mark  LED  lights  for  each 
occurance  or  a  time  pulse  and  stays  lit  for  its  duration. 
The  time  pulse  shows  on  the  chart  paper  as  a  momentary 
increase  in  signal  level  and  is  above  the  local  average 
recorded  signal  strength.  Appendix  C  contains  examples. 
Figure  4.6  is  the  schematic  diagram. 

3 •  Interface  Amplifier 

The  "S"  meter  output  from  the  receiver  and  the  time 
pulse  from  the  timer  are  integrated  into  one  signal  in  the 
interface  amplifier.  This  is  accomplished  by  using  opera¬ 
tional  amplifiers  (LM  148/348).  On  Figure  4.7,  R3  is  used 
as  a  zero  adjust  to  establish  a  zero-signal  reference 
voltage.  hlO  adjusts  the  current  to  the  chart  recorders  so 


the  aeter  stylus  can  be  set  to  its  zero-signal  position.  S5 
is  then  adjusted  for  full  scale  aeter  deflection  at  the 
desired  BF  signal  strength  and  acts  as  a  gain  control. 

4  .  £hapt  fie£.gsflg.£S 

The  fiustrak  recorder  stylus  is  gear  driven  and 
strikes  the  pressure  sensitive  paper  at  regular  intervals  to 
produce  an  apparently  continuous  line.  The  stylus*  position 
depends  upon  the  current  level  through  the  aeter.  The  paper 
speed  is  deterained  by  interchangeable  gear  trains.  Tvo 
styles  of  recorders  are  used  and  both  are  geared  to  allow 
one  rcll  of  paper  to  last  two  weeks.  The  main  functional 
difference  between  the  styles  is  the  aaount  of  current 
required  for  full  scale  deflection. 

5 .  Batter y  Backup 

Figure  4.8  shows  the  power  supply  block  diagraa. 
When  60  HZ,  115  volt  power  is  availaole,  the  12  volt  power 
supply  is  used  to  provide  system  power.  The  relay  contacts 
are  held  open  by  the  60  HZ  power  and  the  battery  chargers 
maintain  full  charge  on  both  batteries.  The  five  diodes  in 
series  <for  both  positive  and  negative  supplies)  are 
reversed  biased  with  no  current  flow.  In  the  event  AC  power 
is  lost,  the  relay  contacts  close  and  power  is  supplied  by 
the  batteries  through  the  three  parallel  diodes.  Three 
diodes  in  parallel  were  used  for  current  handling  and  reli¬ 
ability  reasons.  When  AC  power  is  regained,  the  system 
returns  to  its  original  condition. 

It  was  important  to  maintain  power  to  the  timer 
without  interruption  or  the  time  sequence  became  scrambled. 
This  is  accomplished  by  the  five  series  diodes  for  the  time 
interval  between  losing  AC  power  and  the  relay  contacts 
closing.  The  parallel  diodes  drop  the  battery  voltage  by 
approximately  0.6  volts.  This  was  necessary  because  the 
fully  charged  batteries  produce  12.6  volts  for  several  days. 
While  not  critical  to  the  electronics  portion  of  the  system, 
the  chart  recorder  paper  drive  motors  are  DC  powered  and  the 
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V.  BE SOLIS  AND  CONCLUSIONS 

A.  BESOLTS 

The  system  installation  was  completed  and  began  opera¬ 
tion  on  November  8,  1985.  The  first  set  of  strip  chart 
recordings  cover  the  time  period  of  8-22  November,  1985. 
Appendix  C  contains  samples  from  channel  seven  which  is  at 
150.75  MHz,  vertically  polarized,  and  looking  toward  the 
southeast.  During  the  timeframe  shown,  three  small  earth¬ 
quakes  occured.  The  first  was  a  magnitude  2.3  at  2044  on 
November  11,  the  second  a  magnitude  2.9  at  1411  on 
November  12  and  the  third  a  magnitude  2.04  at  0602  on 
November  13.  No  other  channels  recorded  any  signals  during 
the  same  time  period.  The  magnitude  2.3  was  17  miles  from 
the  antenna  site  with  the  other  two  being  closer.  All  tnree 
were  in  a  southeast  direction  from  the  antennas.  Channel 
ten  (38.45  MHz,  omni)  recorded  two  groups  of  signals  witn 
the  first  occuring  during  the  three  days  preceding  the 
quakes.  The  second  group  was  two  days  after  the  quakes. 


B.  CONCLUSIONS 

While  no  firm  conclusions  can  be  reached  from  such  a 
limited  set  of  chart  recordings,  the  electromagnetic 
activity  on  channels  seven  and  ten  lend  encouragement  to  the 
theory  of  earthquake  generated  RF  energy.  Additional 
recordings  are  required  to  confirm  this  theory  or  show  just 
a  random  coincidence  on  the  first  attempt  at  data 
collection. 

In  order  to  show  a  connection,  recordings  should  be 
continued  for  at  least  a  year.  Not  only  tne  simultaneous 
occurence  but  also  those  where  a  quake  occured  without 
recorded  RF  energy  (and  vice  versa)  would  give  valuable 
information  concerning  reliability  and  false  alarm  rate. 
Any  relationship  between  the  recorded  signal  level  and  the 


magnitude  of  the  quake  or  depth  and  linear  distance  of  the 
epicenter  from  the  antennas  should  also  he  investigated. 
The  present  equipment  will  provide  an  initial  level  of 
confidence  on  which  to  base  a  future  decision  of  whether  or 
not  to  continue  the  research.  If  sufficient  evidence  does 
exist  to  justify  further  time  and  monies,  then  more  sofisti- 
cated  and  sensitive  equipment  will  be  required  to  detail  any 
further  relationships.  The  ideal  system  would  include  auto¬ 
mated  data  collection  similar  to  those  currently  in  use  by 
the  Geological  Survey  Office.  When  sufficiently  understood, 
the  possibility  of  transferring  the  monitoring  to  a 
space-based  sensor  could  finally  be  investigated  with  the 
advantages  mentioned  earlier. 
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